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Abstract 

We analyze rigorously the dynamics of the entanglement between two qubits 
which interact only through collective and local environments. Our approach is 
based on the resonance perturbation theory which assumes a small interaction be- 
tween the qubits and the environments. The main advantage of our approach is 
that the expressions for (i) characteristic time-scales, such as decoherence, disen- 
tanglement, and relaxation, and (ii) observables are not limited by finite times. 
We introduce a new classification of decoherence times based on clustering of the 
reduced density matrix elements. The characteristic dynamical properties such as 
creation and decay of entanglement are examined. We also discuss possible ap- 
plications of our results for superconducting quantum computation and quantum 
measurement technologies. 

1 Introduction 

Entanglement plays a very important role in quantum information processes [1~4] (see 
also references therein). Even if different parts of the quantum system (quantum regis- 
ter) are initially disentangled, entanglement naturally appears in the process of quan- 
tum protocols. This "constructive entanglement" must be preserved during the time 
of quantum information processing. On the other hand, the system generally becomes 
entangled with the environment. This "destructive entanglement" must be minimized 
in order to achieve a needed fidelity of quantum algorithms. The importance of these 
effects calls for the development of rigorous mathematical tools for analyzing the dy- 
namics of entanglement and for controlling the processes of constructive and destruc- 
tive entanglement. Another problem which is closely related to quantum information 
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is quantum measurement. Usually, for a qubit (quantum two-level system), quantum 
measurements operate under the condition Hlo ^ ksT, where T is the temperature, 
CO is the transition frequency, h is the Planck constant, and A;^ is the Boltzmann con- 
stant. This condition is widely used in superconducting quantum computation, when 
T ~ 10 — 20mK and huj/kB ~ 100 — 150mK. In this case, one can use Josephson 
junctions (JJ) and superconducting quantum interference devices (SQUIDs), both as 
qubits [5-9] and as spectrometers [10] measuring a spectrum of noise and other im- 
portant effects induced by the interaction with the environment. Understanding the 
dynamical characteristics of entanglement through the environment on a large time 
interval will help to develop new technologies for measurements not only of spectral 
properties, but also of quantum correlations induced by the environment. 

In this paper, we develop a consistent perturbation theory of quantum dynamics of 
entanglement which is valid for arbitrary times t > 0. This is important in many real 
situations because (i) the characteristic times which usually appear in quantum systems 
with two and more qubits involve different time-scales, ranging from a relatively fast 
decay of entanglement and different reduced density matrix elements (decoherence) to 
possibly quite large relaxation times, and (ii) for not exactly solvable quantum Hamil- 
tonians (describing the energy exchange between the system and the environment) one 
can only use a pcrturbative approach in order to estimate the characteristic dynamical 
parameters of the system. Note, that generally not only are the time-scales for decoher- 
ence and entanglement different, but so are their functional time-dependences. Indeed, 
usually the off-diagonal reduced density matrix elements in the basis of the quantum 
register do not decay to zero for large times, but remain at the level of O(A^), where 
A is a characteristic constant of interaction between a qubit and an environment [11]. 
On the other hand, entanglement has a different functional time dependence, and in 
many cases decays to zero in finite time. Another problem which we analyze in this 
paper is a well-known cut-off procedure which one must introduce for high frequencies 
of the environment in order to have finite expressions for the interaction Hamiltonian 
between the quantum register and the environment. Generally, this artificial cut-off 
frequency enters all expressions in the theory for physical parameters, including decay 
rates and dynamics of observables. At the same time, one does not have this cut-off 
problem in real experimental situations. So, it would be very desirable to develop a 
regular theoretical approach to derive physical expressions which do not include the 
cut-off parameter. We show that our approach allows us to derive these cut-off inde- 
pendent expressions as the main terms of the perturbation theory, which is of O(A^). 
The cut-off terms are included in the corrections of O(A^). At the same time, the 
low-frequency divergencies still remain in the theory, and need additional conditions 
for their removal. 

We describe the characteristic dynamical properties of the simplest quantum regis- 
ter which consists of two not directly interacting qubits (effective spins), which interact 
with local and collective environments. We introduce a classification of the decoherence 
times based on a partition of the reduced density matrix elements in the energy basis 
into clusters. This classification, valid for general A^-level systems coupled to reser- 
voirs, is rather important for dealing with quantum algorithms with large registers. 
Indeed, in this case different orders of "quantumness" decay on different time-scales. 
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The classification of decoherence time-scales which wc suggest will help to separate 
environment-induced effects which are important from the unimportant ones for per- 
forming a specific quantum algorithm. We point out that all the populations (diagonal 
of density matrix) always belong to the same cluster to which is associated the relax- 
ation time. 

We present analytical and numerical results for decay and creation of entangle- 
ment for both solvable (integrable, energy conserving) and unsolvable (non-integrable, 
energy-exchange) models, and explain the relations between them. 

This paper is devoted to a physical and numerical discussion of the dynamical 
resonance theory, and its application to the evolution of entanglement. A detailed 
exposition of the resonance method can be found in [11,12]. As the mathematical 
details leading to certain expressions used in the discussion presented in this paper are 
rather lengthy, we report them separately in [12]. 

2 Model 

We consider two qubits Si and S2, each one coupled to a local reservoir, and both 
together coupled to a collective reservoir. The Hamiltonian of the two qubits is 

i?s = + 52^1, (2.1) 

where Bj = hujj/2 are effective magnetic fields, Uj is the transition frequency, and 5| 
is the Pauli spin operator of qubit j. The eigenvalues of are 

El = Bi -\- B2, E2 = Bi — B2, = —Bi + B2, -E4 = —Bi — B2, (2-2) 

with corresponding eigenstates 

^1 = I + +), ^2 = I + -), ^3 = I - +), ^4 = 1- -), (2.3) 

where 5^|ib) = ±|ib). Each of the three reservoirs consists of free thermal bosons at 
temperature T = 1//3 > 0, with Hamiltonian 

= J2 ^kalkHk, 3 = 1> 2, 3. (2.4) 

k 

The index 3 labels the collective reservoir. The creation and annihilation operators sat- 
isfy [o-j,kj<iji k'] ~ ^jj'^k,k'- The interaction between each qubit and each reservoir has 
two parts: an energy conserving and an energy exchange one. The total Hamiltonian 
is 



H = Hs + Hr, + Hk, + Hr, (2.5) 

+ {XiSf + X2S'i)®cl)3{9) (2.6) 

+ {KiSt + K2SI)C^Mf) (2.7) 

+ fiiSf Mg) + f^2S^ ^ Mg) (2-8) 

+ uiSf»Mf) + ^2SI^Mf)- (2-9) 
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Here, (l)j{g) = ^{a]{g) + aj{g)), with 



= "'^^^ = ^9*k<^j,k- (2-10) 

k k 

The A, K, /X, u are the dimensionless couphng constants. The collective interaction is 
given by (2.6) (energy-exchange) and (2.7) (energy conserving), the local interactions 
are given by (2.8), (2.9). Also, Sj is the spin- flip operator (Pauli matrix) of qubit 
j. In the continuous mode limit, becomes a function g{k), k G M'^. Our approach 
is based on analytic spectral deformation methods [11] and requires some analyticity 
of the form factors f,g. Instead of presenting this condition we will work here with 
examples satisfying the condition. 

(A) Let r > 0, S G 5^ be the spherical coordinates of R^. The form factors h = f,g 
(see (2.6)-(2.9)) are h{r, S) = rPe^^'"/ii(S), with p = -1/2 n, n = 0, 1, . . . and 
m= 1,2. Here, hi is any angular function. 

This family contains the usual physical form factors [13]. We point out that we include 
an ultraviolet cutoff in the interaction in order for the model to be well defined. (The 
minimal mathematical condition for this is that f{k),g{k) be square integrable over 
k G M^.) However, as discussed in point 2. before equation (4.17), our approach 
yields expressions for decay and relaxation rates which, to lowest order in the couplings 
between the qubits and the reservoirs, do not depend on the ultraviolet characteristics 
of the model. 



3 Evolution of qubits: resonance approximation 

We take initial states where the qubits are not entangled with the reservoirs. Let ps 
be an arbitrary initial state of the qubits, and let be the thermal equilibrium state 
of reservoir Rj. Let ps{t) be the reduced density matrix of the two qubits at time t. 
The reduced density matrix elements in the energy basis are 

[ps{t)]mn ■= {^m, Ps{t)^n) 

= IVr,+r,+r3 [ps (8) PR, ® PR2 ® PR3 e-''^/^\0n) {^m\ e'*^/^] , (3.1) 

where we take the trace over all reservoir degrees of freedom. Under the non-interacting 
dynamics (all coupling parameters zero), we have 

= e'''-/^[pmmn, (3.2) 

where Cmn = Em - En- 

In the rest of the paper we use the dimensionless functions and parameters. For 
this we introduce a characteristic frequency, loq, to be defined later, in Section 8, 
and the dimensionless energies, temperature, frequencies and wave vectors of thermal 
excitations, and time by setting 

E'^ = En/{hwo), fk = fk/{n^o), 9k = 9k/{f^o), T' = kBT/{nwo), (3.3) 
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13' = l/T', = Wfe/wo, k' = ck/uo, t' = uot, 

where c is the speed of hght. Below we omit index "prime" in all expressions. 

As the interactions with the reservoirs are turned on (some of kj , \j , jij , Vj nonzero) , 
the dynamics (3.2) undergoes two qualitative changes. 

1. The "Bohr frequencies" 

ee{Ek-Ei :Ek,Eiespec{Hs)} (3.4) 

in the exponent of (3.2) become complex resonance energies, e i— )■ Se, satisfying 
S5£e > 0. If > then the corresponding density matrix elements decay to 
zero (irreversibility). 

2. The matrix elements do not evolve independently any more. To lowest order 
in the couplings, all matrix elements with (m,n) belonging to a fixed energy 
difference E^ — En will evolve in a coupled manner. Thus to a given energy 
difference e, (3.4), we associate the cluster of matrix element indexes 

C{e) = {{k,l) : Ek-Ei = e}. (3.5) 

Both effects are small if the coupling is small, and they can be described by perturbation 
theory of energy differences (3.4). We view the latter as the eigenvalues of the Liouville 
operator 

Ls = Hs^ls-h^Hs, (3.6) 

acting on the doubled Hilbert space "X" (and Tis = C^). The appearance 
of 'complex energies' for open systems is well known to stem from passing from a 
Hamiltonian dynamics to an effective non-Hamiltonian one by tracing out reservoir 
degrees of freedom. The fact that independent clusters arise in the dynamics to lowest 
order in the coupling can be understood hcuristically as follows. The interactions 
change the effective energy of the two qubits, i.e. the basis in which the reduced density 
matrix is diagonal. Thus the eigenbasis of Lg (3.6) is changed. However, to lowest order 
in the perturbation, spectral subspaces with fixed e G spec(Ls) are left invariant and 
stay orthogonal for different unperturbed e. So matrix elements associated to C(e) get 
mixed, but they do not mix with those in C(e'), e 7^ e'. 

Let e be an eigenvalue of Ls of multiplicity mult(e). As the coupling parameters are 
turned on, there are generally many distinct resonance energies bifurcating out of e. We 

(s) 

denote them by £e , where the parameter s distinguishes different resonance energies 
and varies between s = 1 and s = I'ie), where i'{e) is some number not exceeding 
mult(e). We have a perturbation expansion 

4«)=e + 4^) + 0(x^), (3.7) 

where 

X := max{|Kj|, |Aj|, l/iji, jp'jl : j = l,2} (3-8) 
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and where Sj^ = 0{x'^) and 9(5e > 0. The lowest order corrections 5e are the 
eigenvalues of an explicit level shift operator Ag (see [11]), acting on the eigenspace of 
Ls associated to e. There are two bases {tii^} and {rff^} of the eigenspace, satisfying 

Ae# = 5^:h'i\ [KT^:^ = Si'^* {4'\vi''^) = Ss,s'. (3.9) 

(s) ~(s) 

We call the eigenvectors r]e and r]e the 'resonance vectors'. We take interaction pa- 
rameters (/, g and the coupling constants) such that the following condition is satisfied. 

(F) There is complete splitting of resonances under perturbation at second order, i.e., 

(s) 

all the 6e are distinct for fixed e and varying s. 

This condition implies in particular that there arc mult(e) distinct resonance energies 
ei^\ s = 1, . . . , mult(e) bifurcating out of e, so that in the above notation, z/(e) = 
mult(e). Explicit evaluation of ^e''^ shows that condition (F) is satisfied for generic 
values of the interaction parameters (see also (4.12)-(4.16)). 

The following result is obtained from a detailed analysis of a representation of 
the reduced dynamics given in [11], and generalized to the present model with three 
reservoirs. The mathematical details are presented in [12]. 

Result on reduced dynamics. Suppose that Conditions (A) and (F) hold. There 
is a constant xq > such that if x < >cq, then we have for all t > 

[pt]mn= J2 At{m,n;k,l) [po]ki + 0{>^), (3.10) 

{k,l)&CiEm-E„) 

where the remainder term is uniform in t >0, and where the amplitudes Af satisfy the 
Chapman-Kolmogorov equation 

At+r{m,n;k,l) = ^ At{m,n;p,q)Ar{p,q;k,l), (3.11) 

{p,q)eC{Em-En) 

fort,r > 0, with initial condition AQ{m,n; k,l) = S^=k^n=i (Kronecker delta). More- 
over, the amplitudes are given in terms of the resonance vectors and resonance energies 
by 

mult{En-Em) 

At{m,n;k,l)= ^ e'*"^t-i^- (^^ ® ^fc, r?i^;J_^^) (r^g.^^, ^„ ® . 

s=l 

(3.12) 

Remark. The upper bound xq satisfies Xq < const. T, where T is the temperature 
of the reservoirs, [11]. 

We will call the first term on the r.h.s. of (3.10) the resonance approximation 

of the reduced density matrix dynamics. 

Discussion. 1. The result shows that to lowest order in x, and homogeneously in 
time, the reduced density matrix elements evolve in clusters. A cluster is determined 
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by indices in a fixed C(e). Within each cluster the dynamics has the structure of 
a Markov process. Moreover, the transition ampHtudes of this process are given by 
the resonance data. They can be calculated explicitly in concrete models. We have 
therefore a simple approximation of the true dynamics, valid homogeneously in time. 
This is an advantage of the resonance representation. A limitafAon is that this method 
cannot describe the evolution of quantities (averages of observables) which are of the 
order of the square of the coupling parameters, since the error in the approximation is 
of the same order. ^ An illustration of this limitation of the method is the large-time 
behaviour of off-diagonal matrix elements. Gcncrically, all off-diagonals decay to a 
limit having the size O(x^), as i — >■ oo [11]. As soon as a matrix element is of order 
O(x^), the resonance approximation (3.10) cannot resolve its dynamics, since it is of 
the same order as the remainder. 

One of our goals is to describe the evolution of entanglement of the qubits (see 
sections 7, 8). From the above explanations, it is clear that the resonance approxi- 
mation is well suited to describe decay of initial entanglement of qubits (if the initial 
entanglement is much larger than 0(x^)). On the other hand, an initially unentangled 
qubit state will typically become entangled due to the interaction with reservoirs. It is 
expected that the entanglement created may be of the same order as the error in the 
approximation (3.10), and hence the question arises if it is possible to describe this pro- 
cess using the resonance approximation. The answer is positive, as we show numerically 
in section 8: indeed we see that the amount of entanglement created is independent of 
the coupling strength. (The effect of changing x is to shift the time-dependent curve 
of entanglement along the time-axis.) 

2. Cluster classification of density matrix elements. The main dynamics 
partitions the reduced dynamics into independent clusters of jointly evolving matrix 
elements, according to (3.5). Depending on the energy level distribution of the two 
isolated qubits and on the interaction parameters, each cluster has its associated decay 
rate. It is possible that some clusters decay very quickly, while some others stay 
populated for much longer times. The resonance dynamics furnishes us with a very 
concrete recipe telling us which parts of the matrix elements disappear when. This 
reveals a pattern of where in the density matrix quantum properties are lost at which 
time. The same feature holds for an arbitrary A?^-level system coupled to reservoirs, 
[11], and notably for complex systems {N » 1). In particular, this approach may 
prove useful in the analysis of feasibility of quantum algorithms performed on A^-qubit 
registers. We point out that that the diagonal belongs always to a single cluster, namely 
the one associated with e = 0. If the energies of the N-level system are degenerate, 
then some off-diagonal matrix elements belong to the same cluster as the diagonal as 
well. 

3. The sum in (3.10) alone, which is the main term in the expansion, preserves 
the hermiticity but not the positivity of density matrices. In other words, the matrix 

obtained from this sum may have negative eigenvalues. Since by adding O(x^) wc do 
get a true density matrix, the mentioned negativity of eigenvalues can only be of O(x^). 
This can cause complications if one tries to calculate for instance the concurrence by 

^However, by including higher order terms in the perturbation theory, one can refine the resonance 
method and resolve processes of higher order in the coupling. 
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using the main term in (3.10) alone. Indeed, concurrence is not defined in general for 
a 'density matrix' having negative eigenvalues. See also section 8, Numerical Results. 

4. It is well known that the time decay of matrix elements is not exponential 
for all times. For example, for small times it is quadratic in t [13]. How is this 
behaviour compatible with the representation (3.10), (3.12), where only exponential 
time factors e'*^ are present? The answer is that up to an error term of O(x^), the 
"overlap coefHcients" (scalar products in (3.12)) mix the exponentials in such a way as 
to produce the correct time behaviour. 

(s) 

5. Since the coupled system has an equilibrium state, one of the resonances Eq ' is 
always zero [14], we set e''^^ = 0. The condition Qei^^ > for all e, s except e = 0, s = 1 
is equivalent to the entire system (qubits plus reservoirs) converging to its equilibrium 
state for large times. 

As we have remarked above, the decay of matrix elements is not in general expo- 
nential, but we can nevertheless represent it (approximate to order x^) in terms of 
superpositions of exponentials, for all times t > 0. In regimes where the actual dynam- 
ics has exponential decay, the rates coincide with those we obtain from the resonance 
theory (large time dynamics, see Section 5 and also [11,13]). It is therefore reasonable 
to define the thermalization rate by 

7*^ = minSjei'*^ > 

and the cluster decoherence rate associated to C(e), e 0, by 

7^ = min Qei'^ > 0. 

l<s<mult(e) 

The interpretation is that the cluster of matrix elements of the true density matrix 
associated to e 7^ decays to zero, modulo an error term O(x^), at the rate j^^'^, and 
the cluster containing the diagonal approaches its equilibrium (Gibbs) value, modulo 
an error term O(x^), at rate 7*^^. If 7 is any of the above rates, then r = I/7 is the 
corresponding (thermalization, decoherence) time. 

It should be understood that characterizing the dynamcs via the cluster decoherence 
and relaxation times corresponds to a 'coarse graining': matrix elements arc grouped 
into clusters and the dynamics of clusters is effectively given by a decoherence (or 
the relaxation) time. Expression (3.10) gives much more detail, it gives the dynamics 
of each single matrix element. The breakup of the density matrix into individually 
evolving clusters may be advantageous especially in complex systems, where instead of 
two qubits, one deals with A^-qubit registers. 

Remark on the markovian property. In the first point discussed after (3.12), 
we remark that within a cluster, our approximate dynamics of matrix elements has 
the form of a Markov process. In the theory of markovian master equations, one con- 
structs commonly an approximate dynamics given by a markovian quantum dynamical 
semigroup, generated by a so-called Lindblad (or weak coupling) generator [15]. Our 
representation is not in Lindblad form (indeed, it is not even a positive map on density 
matrices). To make the meaning of the markovian property of our dynamics clear. 
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we consider a fixed cluster C and denote the associated pairs of indices by {mk,nk), 
k = 1,...,K. Retaining only the main part in (3.10), and making use of (3.11) we 
obtain for t,s >0 

[Pt+s] mini 



't+s\mKnK 





[Ps\mini 


= Mt) 


. [PslmxnK . 





(3.13) 



where [Ac{t)]mjnj,mini = At{mj,nj;mi,ni), c.f. (3.12). Thus the dynamics of the 
vector having as components the density matrix elements, has the semi-group property 
in the time variable, with generator Gc ■= ^Ac{0), 



[Pt]mini 




[Po\mini 




= ^tGc 




. [Pt\mKnK . 




_ [PolmKnx . 



(3.14) 



This is the meaning of the Markov property of the resonance dynamics. 

While the fact that our resonance approximation is not in the form of the weak 
coupling limit (Lindblad) may represent disadvantages in certain applications, it may 
also allow for a description of effects possibly not visible in a markovian master equa- 
tion approach. Based on results [1,16], one may believe that revival of entanglement 
is a non-markovian effect, in the sense that it is not detectable under the markovian 
master equation dynamics (however, we are not aware of any demonstration of this 
result). Nevertheless, as we show in our numerical analysis below, the resonance ap- 
proximation captures this effect (sec Figure 1). We may attempt to explain this as 
follows. Each cluster is a (indpcndcnt) markov process with its own decay rate, and 
while some clusters may depopulate very quickly, the ones responsible for creating 
revival of entanglement may stay alive for much longer times, hence enabling that pro- 
cess. Clearly, on time-scales larger than the biggest decoherence time of all clusters, the 
matrix is (approximately) diagonal, and typically no revival of entanglement is possible 
any more. 



4 Explicit resonance data 

We consider the Hamiltonian Hg, (2.4), with parameters < Bi < B2 s.t. B2/B1 ^ 2. 
This assumption is a non-degeneracy condition which is not essential for the applica- 
bility of our method (but lightens the exposition). The eigenvalues of Hs are given by 
(2.2) and the spectrum of Ls is {ei, ±62, ±63, ±64, ±65}, with non-negative eigenvalues 

ei = 0, 62 = 2Bi, 63 = 2^2, 64 = 2(^2 - ^1), es = 2(^1 + B2), (4.1) 

having multiplicities mi = 4, m2 = = 2, m.^ = = 1, respectively. According to 
(4.1), the grouping of jointly evolving elements of the density matrix above and on the 
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diagonal is given by^ 





:= C(ei) = 


{(1,1), (2, 2), (3, 3), (4, 4)} 


(4.2) 


C2 


:= C(e2) = 


{(1,3),(2,4)} 


(4.3) 


C3 


:= C(e3) = 


{(1,2), (3, 4)} 


(4.4) 


C4 


:= C(-e4) 


= {(2,3)} 


(4.5) 


C5 


:= C(e5) = 


{(1,4)} 


(4.6) 



There are five clusters of jointly evolving elements (on and above the diagonal). One 
cluster is the diagonal, represented by Ci. For x > and h G L^(M^,d^A;) we define 

c7/,(x) = 47ra;2cothO) / |/i(2x,S)pdS (4.7) 

(spherical coordinates) and for a: = we set 

(7^(0) = 47rlimx^coth(^x) / \h{2x,T,)fdT,. (4.8) 

Furthermore, let 

Y2 = |$5[4K2Kir2-i(Ai + ^2)2^2(^^)_4i^^^^ ^^2^^2) ^^/]l/2|^ ^^ g^ 

F3 = \^[4KlKy-i{\l + i^lfal{Bi)-AiKiK2{Xj + fil)rr[f^\, (4.10) 
(principal value square root with branch cut on negative real axis) where 

r = P.V. / ^d^fc, r'^=^Bj f |5(2S„S)|2dS. (4.11) 

7]R3 \k\ Js2 

The following results are obtained by an explicit calculation of level shift operators. 
Details are presented in [12]. 

Result on decoherence and thermalization rates. The thermalization and 

decoherence rates are given by 





= min{(A2 + ^2^a,(i?,)}+0(x^) 


(4.12) 


72'^^= 


= 5(^1 + /^iK(5i) + 5(^2 + )"i)^9(^2) 






-Y2 + {kI + z/f )cT/(0) + O(x^) 


(4.13) 


73'^'= 


= i(Af + f4)a,{Bi) + i(Ai + fil)ag{B2) 








(4.14) 




= {Xi + ^ll)ag{B,) + {Xl + ^,l)agiB2) 






+ [im - ^^2f + lyf + vf\ (T/(0) + 0(x4) 


(4.15) 


75'''^ 


= {\\-rlx\)a,{B^)-r(\l^ii\)a,(B2) 






+ [(ki + K2f + v\ + (T/(0) + O(X^) 


(4.16) 



^ Since the density matrix is hermitian, it suffices to know the evolution of the elements on and 
above the diagonal. 

10 



Discussion. 1. The thermalization rate depends on energy-exchange parameters \j, 
fij only. This is natural since an energy-conserving dynamics leaves the populations 
constant. If the interaction is purely energy-exchanging {kj = Uj = 0), then all the rates 
depend symmetrically on the local and collective interactions, through X'j+fij. However, 
for purely energy-conserving interactions (Xj = fij = 0) the rates are not symmetrical in 
the local and collective terms. (E.g. 74'^'^ depends only on local interaction if ki = K2-) 
The terms Y2, I3 are complicated nonlinear combinations of exchange and conserving 
terms. This shows that effect of the energy exchange and conserving interactions are 
correlated. 

2. We see from (4.7), (4.8) that the leading orders of the rates (4.12)-(4.16) do not 
depend on an ultraviolet features of the form factors f,g. (However, crjp(O) depends 
on the infrared behaviour.) The coupling constants, e.g. in (4.12) multiply ag{Bj), 
i.e., the rates involve quantities like (see (4.7)) 

TrAf / coth(/3|A:|/2) \g{\k\,i:)\^ 5^'^\\k\ - 2Bj) d^k. (4.17) 

The one-dimensional Dirac delta function appears due to energy conservation of pro- 
cesses of order x^, and 2Bj is (one of) the Bohr frequencies of a qubit. Thus energy 
conservation chooses the evaluation of the form factors at finite momenta and thus an 
ultraviolet cutoff is not visible in these terms. Nevertheless, we do not know how to 
control the error terms O(x^) in (4.12)-(4.16) homogeneously in the cutoff. 

3. The case of a single qubit interacting with a thermal bose gas has been extensively 
studied, and decoherence and thermalization rates for the spin-boson system have been 
found using different techniques, [17-19]. We recover the spin-boson model by setting 
all our couphngs in (2.5)-(2.9) to zero, except for Xi = ki = A, and setting f = g. In 
this case, the spectral density J{oj) of the reservoir is linked to our quantity (4.7) by 



coth(;Sa;/2)' 



The relaxation rate is 



f h = Kx'^coth{pB)J{2B), 



where 2B is the transition frequency of qubit (in units where H = 1), see (2.1). The 
decoherence rate is given by 

where (T/i(0) is the limit as a; ^ of coth(/3a;) J(2a;). These rates obtained with our 
resonance method agree with those obtained in [17-19] by the standard Bloch-Redfield 
approximation. 

Remark on the limitations of the resonance approximation. As mentioned in 
section 3, the dynamics (3.10) can only resolve the evolution of quantities larger than 
O(x^). For instance, assume that in an initial state of the two qubits, all off-diagonal 
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density matrix elements are of the order of unity (relative to x). As time increases, 

J. dec rt. 

the off-diagonal matrix elements decrease, and for times t satisfying e '^^ < 0{>c ), 
the off-diagonal cluster Cj is of the same size 0{x'^) as the error in (3.10). Hence the 
evolution of this cluster can be followed accurately by the resonance approximation for 



times t < ln(x-2)/'ydec 

oc ^^^'^^ temperature. Here, T, x (and other 

parameters) are dimensionless. To describe the cluster in question for larger times, 
one has to push the perturbation theory to higher order in x. It is now clear that 
if a cluster is initially not populated, the resonance approximation does not give any 
information about the evolution of this cluster, other than saying that its elements will 
be O(x^) for all times. 

Below we investigate analytically decay of entanglement (section 6) and numeri- 
cally creation of entanglement (section 8). For the same reasons as just outlined, an 
analytical study of entanglement decay is possible if the initial entanglement is large 
compared to O(x^). However, the study of creation of entanglement is more subtle 
from this point of view, since one must detect the emergence of entanglement, presum- 
ably of order O(x^) only, starting from zero entanglement. We show in our numerical 
analysis that entanglement of size 0.3 is created independently of the value of x (rang- 
ing from 0.01 to 1). We are thus sure that the resonance approximation does detect 
creation of entanglement, even if it may be of the same order of magnitude as the 
couplings. Whether this is correct for other quantities than entanglement is not clear, 
and so far, only numerical investigations seem to be able to give an answer. As an ex- 
ample where things can go wrong with the resonance approximation we mention that 
for small times, the approximate density matrix has negative eigenvalues. This makes 
the notion of concurrence of the approximate density matrix ill-defined for small times. 



5 Comparison between exact solution and resonance ap- 
proximation: explicitly solvable model 

We consider the system with Hamiltonian (2.5)-(2.10) and Ai = A2 = 0, f^i = = 0, 
Ki = K2 = K and vi = 1/2 = 1^- This energy-conserving model can be solved explicitly 
[11,13] and has the exact solution 



where 



and 



[Pt]mn = [Po] 



-lt{Em-En) ^iK^amnSit) g- [«;^6mn+l^^Cmn]r(t) 
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S{t) 

m 



= [ \f{k)fcothiP\k\/2) 



sin^(|fc|V2) 3 



(5.1) 



4 4 8 

4 8 4 

4 8 4 

8 4 4 



(5.2) 
(5.3) 
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On the other hand, the main contribution (the sum) in (3.10) yields the resonance 
approximation to the true dynamics, given by 



[Pt]mm 
[Pt]ln 
[ft] 14 
[Pt] 23 
[Pt]m4 



In 



n 



2,3 



[po]mm m = 1,2,3,4 
g-it(Ei-i?4) e-*(^«'+2-'W(o)[po]i4 

^-itiE,-E,) e-2*-^-/(0)[po]23 

^-it{Em-E,) g2it«2, ^_i(«2+,2)^^(o) [p^]^^ m = 2, 3 



(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 



The dotted equaHty sign = signifies that the left side equals the right side modulo an 
error term 0{k'^ + z/^), homogeneously in t > 0.^ Clearly the decohercncc function 
r(t) and the phase S{t) arc nonlinear in t and depend on the ultraviolet behaviour of 
/. On the other hand, our resonance theory approach yields a representation of the 
dynamics in terms of a superposition of exponentially decaying factors. From (5.1) and 
(5.4)-(5.8) we see that the resonance approximation is obtained from the exact solution 
by making the replacements 



S{t) ^ \rt, 



(5.9) 
(5.10) 



We emphasize again that, according to (3.10), the difference between the exact solu- 
tion and the one given by the resonance approximation is of the order 0{k^ + u'^), 
homogeneously in time, and where O^k^ + v'^) depends on the ultraviolet behaviour of 
the couplings. This shows in particular that up to errors of 0{k^ + v'^), the dynamics 
of density matrix elements is simply given by a phase change and a possibly decay- 
ing exponential factor, both linear in time and entirely determined by r and cr/(0). 
Of course, the advantage of the resonance approximation is that even for not exactly 
solvable models, we can approximate the true (unknown) dynamics by an explicitly 
calculable superposition of exponentials with exponents linear in time, according to 
(3.10). Let us finally mention that one easily sees that 



lim S{t)/t 



r/2 and 



lim T{t)/t 



a/(0)/4. 



so (5.9) and (5.10) may indicate that the resonance approximation is closer to the true 
dynamics for large times - but nevertheless, our analysis proves that the two are close 
together {0{k^ + v'^)) homogeneously in i > 0. 



6 Disentanglement 

In this section we apply the resonance method to obtain estimates on survival and 
death of entanglement under the full dynamics (2.5)-(2.9) and for an initial state of the 
form ps <S> PRi <S> PRj "X" Pr^ , where ps has nonzero entanglement and the reservoir initial 

^To arrive at (5.4)-(5.8) one calculates the At in (3.10) explicitly, to second order in k and v. The 
details are given in [12]. 
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states are thermal, at fixed temperature T = > 0. Let p be the density matrix of 
two qubits 1/2. The concurrence [20,21] is defined by 

C(p) = max{0, [v^ + V^ + -v/z^]}' (6-1) 

where i^i > ^2 > t's > > are the eigenvalues of the matrix 

i{p) = p{sy sy)p*{sy o s^). (6.2) 

Here, p* is obtained from p by representing the latter in the energy basis and then 

taking the elementwise complex conjugate, and is the Pauli matrix S'^ = ^ ^ . 

The concurrence is related in a monotone way to the entanglement of formation, and 
(6.1) takes values in the interval [0,1]. If C(p) = then the state p is separable, 
meaning that p can be written as a mixture of pure product states. If C{p) = 1 we call 
p maximally entangled. 

Let po be an initial state of S. The smallest number to > s.t. C{pt) = 
for all t > to is called the disentanglement time (also 'entanglement sudden death 
time', [1,22]). If C{pt) > for all t > then we set to = oo- The disentanglement time 
depends on the initial state. Consider the family of pure initial states of S given by 



Po 



with 



ai 



+ +> + 



02 



VI«lP + l«2|2' VI«lP + l«2|2 

where ai,a2 G C are arbitrary (not both zero). The initial concurrence is 

|Kaia2| 



C{po 



joij + |a2| 



which covers all values between zero (e.g. ai = 0) to one (e.g. ai = 02 G 
to (3.10), the density matrix of S at time t > is given by 



Pt 



pi a 

P2 

P3 

a* p4 



+ 0(x2). 



According 



(6.3) 



with remainder uniform in t, and where pj = Pjit) and a = a{t) are given by the 
main term on the r.h.s. of (3.10). The initial conditions are pi(0) = , P2(0) = 

^,3(0) = 0, p4(0) = and a(0) = We set 

p:=pi(0) e [0.1] (6.4) 



and note that ^4(0) = 1 — p and |a(0)| = ■\/p{\ — p)- In terms of p, the initial concur- 
rence is C{pq) = 2a/p(1 — p). Let us set 



^2 := (A? + pl)ag{B{), S3 := {Xj + pl)ag{B2), 



(6.5) 
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S5 ■= S2 + S3+ [{ki + K2f + fl + yl] (7/(0). (6.6) 

(5+ := max{ (52,(53}, (5_ := inin{(52, ^a}. (6.7) 

An analysis of the concurrence of (6.3), where the Pj{t) and a(t) evolve according to 
(3.10) yields the following bounds on disentanglement time. 

Result on disentanglement time. Take p 7^ 0, 1 and suppose that 82^8^ > 0. 

There is a constant xq > (independent of p) such that we have: 

A. (Upper bound.) There is a constant Ca > (independent ofp, k) s.t. C{pt) = 
for all t > tA, where 



tA '■= max 



In 



Ca 



VpC^-p) 



'52 + ^3 



In 



Ca 



p{\ -p) 



Ca 



' (52 + (53 I 



(6.8) 



B. (Lower bound.) There is a constant Cb > (independent of p, k) s.t. C{pt) > 
for all t <tB, where 



(6.9) 



Bounds (6.8) and (6.9) are obtained by a detailed analysis of (6.1), with p replaced 
by pt, (6.3). This analysis is quite straightforward but rather lengthy. Details are 
presented in [12]. 

Discussion. 1. The result gives disentanglement bounds for the true dynamics of the 
qubits for interactions which are not integrable. 

2. The disentanglement time is finite. This follows from (52,(53 > (which in turn 
implies that the total system approaches equilibrium as t — )■ 00). If the system does 
not thermalize then it can happen that entanglement stays nonzero for all times (it 
may decay or even stay constant) [1,23]. 

3. The rates 6 are of order x^. Both tA and tB increase with decreasing coupling 
strength. 

4. Bounds (6.8) and (6.9) are not optimal. The disentanglement time bound (6.8) 
depends on both kinds of couplings. The contribution of each interaction decreases 
tA (the bigger the noise the quicker entanglement dies). The bound on entanglement 
survival time (6.9) does not depend on the energy-conserving couplings. 



7 Entanglement creation 

Consider an initial condition ps PB.2 PR3 1 where ps is the initial state of the 
two qubits, and where the reservoir initial states are thermal, at fixed temperature 
T=l/P> 0. 

Suppose that the qubits are not coupled to the collective reservoir R3, but only to 
the local ones, via energy conserving and exchange interactions (local dynamics). It is 
not difficult to see that then, if ps has zero concurrence, its concurrence will remain 
zero for all times. This is so since the dynamics factorizes into parts for Si + Ri and 
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§2 +1^21 and acting upon an unentangled initial state does not change entanglement. 
In contrast, for certain entangled initial states ps, one observes death and revival of 
entanglement [16]: the initial concurrence of the qubits decreases to zero and may 
stay zero for a certain while, but it then grows again to a maximum (lower than the 
initial concurrence) and decreasing to zero again, and so on. The interpretation is that 
concurrence is shifted from the qubits into the (initially unentangled) reservoirs, and 
if the latter are not Markovian, concurrence is shifted back to the qubits (with some 
loss) . 

Suppose now that the two qubits arc coupled only to the collective reservoir, and 
not to the local ones. Braun [24] has considered the explicitly solvable model (energy- 
conserving interaction), as presented in Section 5 with k = 1, u = 0.^ Using the exact 
solution (5.1), Braun calculates the smallest eigenvalue of the partial transpose of the 
density matrix of the two qubits, with S and T considered as non-negative parameters. 
For the initial product state where qubits 1 and 2 are in the states ;^(|+) ~ |~)) a^id 
^(1+) + |-)) respectively, i.e.. 



it is shown that for small values of T (less than 2, roughly), the negativity of the smallest 
eigenvalue of the partial transpose oscillates between zero and -0.5 for S increasing from 
zero. As T takes values larger than about 3, the smallest eigenvalue is zero (regardless 
of the value of S). According to the Peres-Horodecki criterion [25,26], the qubits are 
entangled exactly when the smallest eigenvalue is strictly below zero. Therefore, taking 
into account (5.2) and (5.3), Braun's work [24] shows that for small times (T small) the 
collective environment (with energy-conserving interaction) induces first creation, then 
death and revival of entanglement in the initially unentangled state (7.1), and that for 
large times (F large), entanglement disappears. 

Resonance approximation. The main term of the r.h.s. of (3.10) can be 
calculated explicitly, and we give in Appendix A the concrete expressions. How does 
concurrence evolve under this approximate evolution of the density matrix? 

(1) Purely energy-exchange coupling. In this situation we have k = v = 0. The 
explicit expressions (Appendix A) show that the density matrix elements [pf]mri in the 
resonance approximation depend on A (collective) and /x (local) through the symmet- 
ric combination -|- //^ only. It follows that the dominant dynamics (3.10) (the true 
dynamics modulo an error term O(x^) homogeneously in t > 0) is the same if we 
take purely collective dynamics (^u = 0) or purely local dynamics (A = 0). In particu- 
lar, creation of entanglement under purely collective and purely local energy -exchange 

"'in fact, Brown uses this model and sets the Hamiltonian of the qubits equal to zero. This has no 
influence on the evolution of concurrence, since the free dynamics of the qubits can be factored out 
of the total dynamics (energy-conserving interaction), and a dynamics of Si and S2 which is a prouct 
does not change the concurrence. 




1 



1 

1 
-1 
-1 



1 

1 
-1 
-1 



-1 

-1 
1 
1 



-1 

-1 
1 
1 



(7.1) 
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dynamics is the same, modulo 0{>c'^). For instance, for the initial state (7.1), collec- 
tive energy-exchange couplings can create entanglement of at most O(x^), since local 
energy-exchange couplings do not create any entanglement in this initial state. 

(2) Purely energy- conserving coupling. In this situation we have X = fi = 0. The 
evolution of the density matrix elements is not symmetric as a function of the cou- 
pling constants k (collective) and u (local). One may be tempted to conjecture that 
concurrence is independent of the local coupling parameter ly, since it is so in absence 
of collective coupling {k = 0). However, for k, ^ 0, concurrence depends on (see 
numerical results below). We can understand this as follows. Even if the initial state 
is unentangled, the collective coupling creates quickly a little bit of entanglement and 
therefore the local environment does not see a product state any more, and starts 
processes of creation, death and revival of entanglement. 

(3) Full coupling. In this case all of k, A, fi, v do not vanish. Matrix elements evolve 
as complicated functions of these parameters, showing that the effects of different in- 
teractions are correlated. 



8 Numerical Results 

In the following, we ask whether the resonance approximation is sufficient to detect 
creation of entanglement. To this end, we take the initial condition (7.1) (zero concur- 
rence) and study numerically its evolution under the approximate resonance evolution 
(Appendices A, B), and calculate concurrence as a function of time. Let us first con- 
sider the case of purely energy conserving collective interaction, namely A = // = = 
and only k / 0. Our simulations (Figure la) show that, a concurrence of value ap- 




Figure 1: Energy conserving collective interaction \ = n = v = Q. a) Concurrence a function of 
time for different k values as indicated in the legend, b) The same as a) but in the renormalized 
time K^t. 



17 



proximately 0.3 is created, independently of the value of k (ranging from 0.01 to 1). 
It is clear from the graphs that the effect of varying k consists only in a time shift. 
This shift of time is particularly accurate, as can be seen in Fig. lb, where the three 
curves drawn in a) collapse to a single curve under the time rescaling t — t- K^t. In 
particular, the maximum concurrence is taken at times to ~ 0.5k~^. We also point out 
that the revived concurrence has very small amplitude (approximately 15 times smaller 
than the maximum concurrence) and takes its maximum at ti ~ 2.1k~^. Even though 
the amplitude of the revived concurrence is small as compared to k^, the graphs show 
that it is independent of k, and hence our resonance dynamics does reveal concurrence 
revival. 




Figure 2: Energy conserving collective and local interaction \ = ji = Q. a) Concurrence a 
function of time for fixed collective interaction k = O.Ol and different local interaction v as 
indicated in the legend, b) Variation of the maximum of concurrence as a function of the local 
interaction strength v for different collective interaction strengths k as indicated in the legend. 

When switching on the local energy conserving coupling, 7^ 0, we see in Fig. 2a, 
that the maximum of concurrence decreases with increasing v. Therefore, the effect 
of a local coupling is to reduce the entanglement. It is also interesting to study the 
dependence of the maximal value of the function of the energy- 

conserving interaction parameters. This is done in Fig. 2b, where Cmax is plotted as a 
function of the local interaction for different fixed collective couplings k. The graphs 
show that as the local coupling v is increased to the value of the collective coupling 

Cmax becomes zero. This means that if the local coupling exceeds the collective 
one, then there is no creation of concurrence. We may interpret this as a competition 
between the concurrence-reducing tendency of the local coupling (apart from very small 
revival effects) and the concurrence-creating tendency of the collective coupling (for not 
too long times). If the local coupling exceeds the collective one, then concurrence is 
prevented from building up. 

Looking at Fig. 2, it is clear that the effect of the local coupling is not only to 
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Figure 3: Energy conserving collective and local interaction X = ^ = 0. Rescaled concurrence 
C{p) / Cmax as function of time for fixed local interaction v = 0.005 and different collective 
interaction k > ly (as indicated in the legend) as a function of the rescaled time (k^ + j/^)t. 




Figure 4: Energy-exchanging collective and local interactions X — fj, 0. Concurrence C{p) as 
function of time for fixed energy-conserving collective interactions k = 0.02, ly = and different 
energy-exchanging couplings A as indicated in the legend. Here we used Bi = I, B2 = 1.25, 
and (3 — 1. 

decrease concurrence but also to induce a shift of time, similarly to the effect of the 
collective coupHng k. Indeed, taking as a variable the rescaled concurrence C^p)/Crnax^ 
one can see that the approximate scaling (k^ -|- z^^)t is at work, see Fig. 3. We conclude 
that both local and collective energy conserving interactions produce a cooperative time 
shift of the entanglement creation, hut only the local interaction can destroy entangle- 
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ment creation. There is no entanglement creation for v > k. 

Let us now consider an additional energy exchange coupling A, ^ / 0. Since these 
parameters appear in the resonance dynamics only in the combination + /i^, see 
Appendix A, we set without loosing generality A = We plot in Fig. 4 the time 
evolution of the concurrence, at fixed energy-conserving couplings k = 0.02 and = 0, 
for different values of the energy exchange coupling A. In this case we have chosen 
Bi = 1 which corresponds to ujq = uJi/2, where loi is a transition frequency of the 
first qubit. We also used the conditions: ag(Bi) = rg(Bi) = 1, which lead to the 
renormalization of the interaction constants. The relations between (Tg{B2) and ag{Bi), 
and rg{B2) and rg{Bi) are discussed in Appendix B. 

Figure 4 shows that the effect of the energy exchanging coupling is to shift slightly 
the time where concurrence is maximal and, at the same time, to decrease the amplitude 
of concurrence for each fixed time. This feature is analogous to the effect of local 
energy-conserving interactions, as discussed above. Unfortunately, it is quite difficult 
in this case to extract the threshold values of A at which the creation of concurrence 
is prevented for all times. The difficulty comes from the fact that for larger values 
of A, the concurrence is very small and the negative eigenvalues on order 0{>c^) do 
not allow a reliable calculation. This picture docs not change much if a local energy- 
conserving interaction v < k \s added. In Fig. 5, we show respectively, the time shift 
of the maximal concurrence Ai = tmaxW — tmax{^ = 0) as a function of the energy- 
exchanging coupling A (a) and the behavior of the maximal concurrence as a function 
of the same parameter A for two different values of the local coupling i'. Is appears 
evident that the role played by the energy-exchange coupling is very similar to that 
played by the local energy-conserving one. 

Let us comment about concurrence revival. The effect of a collective energy- 
conserving coupling consists of creating entanglement, destroying it and creating it 
again but with a smaller amplitude. Generally speaking, an energy-exchanging cou- 
pling, if extremely small, does not change this picture. Nevertheless, it is important 
to stress that the damping effect the energy-exchange coupling has on the concurrence 
amplitude is stronger on the revived concurrence than on the initially created one. 
This is shown in Fig. 6, where the renormalized concurrence C{p) /Cmax is plotted for 
different A values. For these parameter values, only a very small coupling A < 0.001 
will allow revival of concurrence. 

In the calcTilation of concurrence, the sqTiarc roots of the eigenvalues of the matrix 
^{p) (6.2) should be taken. As explained before, the non positivity, to order O(x^) of 
the density matrix p reflects on the non positivity of the eigenvalues of the matrix ^{p). 
When this happens (fj < 0) we simply put i^j = in the numerical calculations. This 
produces an approximate (order 0(x^)) concurrence which produces spurious effects, 
especially for small time, when concurrence is small. These effects are particularly evi- 
dent in Fig. 6, for small time, where artificial oscillations occur, instead of an expected 
smooth behavior. In contrast to this behaviour, the revival of entanglement as revealed 
in Figure 6 varies smoothly in A, indicating that this effect is not created due to the 
approximation. 
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Figure 5: Energy-exchanging collective and local interaction \ = jj, ^Q. a) Time shift induced 
by energy-exchanging coupling, for the same energy conserving collective coupling n = 0.02 and 
different local couplings v as indicated in the legend, b) Decay of the maximal concurrence as 
a function of A, for the same cases as (a). Magnetic fields and temperature is the same as in 
Fig. 4. 




Figure 6: Energy-exchanging collective and local interaction A = /i 7^ 0. Rescaled concurrence 
C{p) / Cmax vs time t, for different A values. Here,K = 0.02 and — 0. Magnetic fields and 
temperature is the same as in Fig. 4. 

9 Conclusion 

We consider a system of two qubits interacting with local and collective thermal quan- 
tum reservoirs. Each qubit is coupled to its local reservoir by two channels, an energy- 
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conserving and an energy-exchange one. The qubits are cohectively coupled to a third 
reservoir, again through two channels. This is thus a versatile model, describing local 
and collective, energy-conserving and energy-exchange processes. 

We present an approximate dynamics which describes the evolution of the reduced 
density matrix for all times t > 0, modulo an error term O(x^), where k is the typi- 
cal coupling strength between a single qubit and a single reservoir. The error term is 
controlled rigorously and for all times. The approximate dynamics is markovian and 
shows that different parts of the reduced density matrix evolve together, but indepen- 
dently from other parts. This partitioning of the density matrix into clusters induces a 
classification of decoherence times - the time-scales during which a given cluster stays 
populated. We obtain explicitly the decoherence and relaxation times and show that 
their leading expressions (lowest nontrivial order in x) is independent of the ultravi- 
olet behaviour of the system, and in particular, independent of any ultraviolet cutoff, 
artificially needed to make the models mathematically well defined. 

We obtain analytical estimates on entanglement death and entanglement survival 
times for a class of initially entangled qubit states, evolving under the full, not explic- 
itly solvable dynamics. We investigate numerically the phenomenon of entanglement 
creation and show that the approximate dynamics, even though it is markovian, does 
reveal creation, sudden death and revival of entanglement. We encounter in the nu- 
merical study a disadvantage of the approximation, namely that it is not positivity 
preserving, meaning that for small times, the approximate density matrix has slightly 
negative eigenvalues. 

The above-mentioned cluster-partitioning of the density matrix is valid for general 
A''-level systems coupled to reservoirs. We think this clustering will play a useful and 
important role in the analysis of quantum algorithms. Indeed, it allows one to separate 
"significant" from "insignificant" quantum effects, especially when dealing with large 
quantum registers for performing quantum algorithms. Depending on the algorithm, 
fast decay of some blocks of the reduced density matrix elements can still be tolerable 
for performing the algorithm with high fidelity. 

We point out a further possible application of our method to novel quantum measur- 
ing technologies based on superconducting qubits. Using two superconducting qubits 
as measuring devices together with the scheme considered in this paper will allow one 
to extract not only the special density of noise, but also possible quantum correlations 
imposed by the environment. Modern methods of quantum state tomography will allow 
to resolve these issues. 

A Dynamics in resonance approximation 

We take < Bi < B2, B2/B1 / 2, and « min{25i,2(B2 - Bi),2\B2 - 2Bi\}. 
These conditions guarantee that the resonances do not overlap, see also [11]. In the 
sequel, = means equality modulo an error term O(x^) which is homogeneous in t > 0. 
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The main contribution of the dynamics in (3.10) is given as follows. 
1 1 
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Of course, the populations do not depend on any energy- conserving parameter. The 
cluster of matrix elements {(3, 1), (4,2)} evolves as 
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^^^^i^{H3i +y-[po]42}. (A.ll) 
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Here, 



1/2 



e'^>^ =A+ -B{1 + es) - (-1)'^- [B\l + ea)^ + AC{B{e2 - 1) + C)J , (A.12) 
where 



A = i{x■' + ^J,^)^ag{B,) + i{^ + l^^)af{0)-{\^ + ^lyg{B^) (A.13) 

B = + ^^)a-{B2) (A.14) 

C = -2kV/ (A.15) 
A + C — e^''^ 

y± = 1 + (A; = 1 for y+, k = 2 for (A.16) 

62-0 



and 



C7g(x) = Attx'^ coth{Px) / |£?(2x,S)pdS 

^ 2..^-^/^J,(2.),E)PdE 

r,(x) = iP.V./j,^5,n2|5(|n|,S)pcoth(^|n|/2)^dndS (A.17) 

r/ = P.V. / ^d^A; 

Jrs \k\ 

(7/(0) = 47rlimx2coth(;9x) / |/(2ar), S)pdS. 
The cluster of matrix elements {(2, 1), (4,3)} evolves as 

[pt]21 = e 2^2-__^-^|[po]21+2/+[po]43j 

+^ ^""^ 1 + ei{y' )2 + ^-[Po]^^}, (A.18) 

[pt]43 = e*^2i2__^^t_|[pQ]2^+y^[pjj]^3| 

i:p^^{[Po]2i + 2/'_[po]43}. (A.19) 

(k) (k) 

Here, £2^2 same as £2^1 ' but with all indexes labeling qubits 1 and 2 interchanged 

(ei O 62, -Bi o in all coefficients involved in £3^^ above). Also, y'^ is obtained 
from ?/± by the same switch of labels. Finally, 

[Pth2 = e'*^2(si-B2)[pQ]32 (A.20) 

[ptUi = e'*^2(«i+«2)[po]4i (A.21) 
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with 



+(A2 + A.2)[r,(Si)-r,(S2)] 
-(A2 + /x2)[r,(Si)+r,(B2)]. 



B Reduction to independent parameters 

The equations above contain four independent coupUng constants A, k describing 
the energy-conserving and the energy exchanging (local and collective) interaction, and 
eight different functions of the form factors / and g : ag{Bi), rg{Bi), ag{Bi), i = 1,2, 
af{0), rf (A.17). 

These functions are not independent. First of all it is easy to see that the following 
relation holds: 

moreover, choosing for instance a form factor g{2x, S) oc \/2x one has: 

(7^(52) _ ^i?2^'coth/3i?2 ^3 2) 



(Tg(5i) \BiJ coih. 15 Bi 

Integrals in du in Eq. (A.17) converge only when adding a cut-off Uc- It is easy to show 
that, when Uc ^ oo one has: 

lim ^ = 1, (B.3) 



uc-)-oo rg{Bi) 



and we can assume rg{Bi) ~ rg{B2). So, we end up with four independent divergent 
integrals, ag{B\), rg{B\), cr/(0), r/, in terms of which we can write explicitly the decay 
rEtijGS * 

ai = {\^ + ^?)ag{Bx) 



(B.4) 



and the Lamb shifts. 



a, = (A^ + /x>.(5i)(t)'iS|t 

as = K^cr/(0) 

0:4 = l'^(T/(0), 

/3i = {X^ + iJ,^)rg{Br) 

/32 = (X'' + n'')rgiB2) Pi (B.5) 

P3 = -K^rf. 
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Suppose now that both Lamb shifts, and decay constants are experimentally measur- 
able quantities, and also assume (due to symmetry) that \ = jx. Interaction constants 
can be renormalized in order to give directly decay constants and Lamb shifts: 

ai = 2A2 

2 ( B2\^ coth/3B2 



«2 = 2V(f , 



as = 



0:4 = 



(B.6) 



h = 2A2 

132 = /5i 

A, K,D are the values chosen for simulations. 
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